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The Schwinger-Dyson equation for a scalar propagator is solved in Minkowski space with 
the help of an integral spectral representation, both for spacelike and timelike momenta. 
The equation is re-written into a form suitable for numerical solution by iterations. This 
procedure is described for a simple unphysical Lagrangian with a cubic interaction, with 
future extensions to more realistic theories in mind. 

1. INTRODUCTION 

The infinite tower of integral Schwinger-Dyson equations (SDEs) links n-point Green's 
functions {n = 2, 3 . . .) of a quantum-field theory. Their exact solution would provide 
complete information on the physics of the theory, including its non-perturbative regime. 
In practice, the system of SDEs has to be truncated and closed by making assumptions 
about the driving term, which contains more complicated Green's functions determined 
by equations that were thrown away. Most often only the simplest equation for the 2- 
point Green's function(s) (one particle propagator) is retained, its kernel contains n-point 
(n = 3, . . .) vertex the full form of which is unknown. Hence some physical ansatz has 
to be used and then the equations can be solved (usually numerically). One hopes that 
the solution provides some useful information on the behavior of the theory, in particular 
in the non-perturbative region. It is certainly interesting to compare the solution with 
results of alternative non-perturbative techniques 

In most papers on the solution of the SDK, the Wick rotation from Minkowski to 
Euclidean space is employed in order to escape singularities of the kernel inherent to 
physical Green functions. We instead attempt to find the solution directly in Minkowski 
space by making use of the spectral representation of the Green's functions based on 
their known (or assumed) analytical properties. In this contribution we will for simplicity 
illustrate this method of solution with a simple example of self-interacting scalar fields with 
super- renormalizable unphysical g^^{x) coupling. Possibilities of extending the technique 
to more realistic theories are briefly discussed. 
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2. FORMALISM 

Our approach is a straightforward extension of the spectral decomposition method of 
Kusaka et al. [Q, developed for the solution of the Bethe-Salpeter equation for scalar 
bound states. The Green's functions - one scalar particle propagator in our case - are 
written as spectral integrals over some weight functions. Then we put these parameter- 
izations and an expression for the vertex function into the SDE, combine denominators 
with the help of the usual Feynman parameterization, integrate over the loop momen- 
tum/momenta and obtain the real integral equations for the weight functions depending 
on spectral variables. These equations are free of singularities and can be solved by 
iterations. 

The generic spectral decomposition of the dressed renormalized scalar propagator reads: 

G{p^) = [ da^r^^ , a{a) = 5{a - m^) + p{a) , (1) 

J — a — ie 

where a{a) is a positive spectral function. Here, we assume that a has a singular contri- 
bution due to the propagation of physical particle with the mass m and a regular positive 
smooth part p{a) which starts at the two-particle threshold. This assumption means that 
the particle spectrum of the system is essentially perturbative, e.g., there is no confine- 
ment and also the contribution of possible bound states below ath = 2m is neglected (the 
latter appears as corrections to the vertex from the poles in the 4-point Green's function). 

The form of the SDE for an one particle propagator depends on the form of the inter- 
action. For our toy model with Cmt = —g^'^{x) and introducing the renormalized (by 
on- mass shell subtraction) self-energy n/j(p^) we get 



G-\p^) = z{p^ - ml) - n(p2) = / _ ^2 _ n^^^2) . n^^^2) ^ ^n^(m2) = 02) 

n(p') = ^7l4 /Ar(p,g)G((p-g)^)G(g2), (3) 



(2vr)4 . 

where mo is a bare mass, the residue of G{p'^) at the pole p^ = m^ equals 1; 5* = 18 is a 
combinatorial factor, Z is the renormalization function of the $(x) field, which is finite 
for the cubic interaction. The adopted renormalization procedure suggests the following 
spectral decomposition for IljiijP'): 

^r{p') = ip' - [ da /-^^^ . . (4) 

J p'^ — a — le 

We will mostly work in the bare vertex approximation r(p, g) = 1. From the imaginary 
part of G = Go + GqUrG and from the dispersion relations of Eqs. (|l|) and (|D follows: 

I ^ I ^ , I 2 \TD r A P7r(^)p(a) +P^(«)P(^) .^x 

p[uj) = Ptt[uj) + [m — uj) P da . (5j 

J4m2 a — LU 

From the SDE (^ we get the second relation between and p: 
Sg' 1 



P^(^) 



(27r)4 {uj-m^f 



•y/l — Am^ /uj + / daT{a,m^ ,uj)p{a) 



+ 



j J daida2T{ai,a2,uj)p{ai)p{a2) 



(6) 



3 



where T{ai, a2, uj) is a purely kinematical function which also determines the integration 
bounds. The set of equations is solved by iterations. 

The method works in essentially the same way also for more complicated theories. Of 
course, if there is more than one type of particles, one gets a set of coupled integral 
equations. In more complicated theories it might not be always possible to integrate out 
all Feynman parameters. Then kinematical functions analogous to T{ai, 02,00) in Eq. 
(P) are expressed in terms of the remaining integrations. That might slow down and 
complicate the solution, in particular one has to determine numerically the bounds for 
integrals over the spectral parameters a. 

3. RESULTS 

The iterations work very well for small A = S g"^ / {{2'kY im?) , but in that region the 
propagator is not far from the free one. With increasing A the rate of convergence is 
slowing down rather fast, then for a certain value \crit ~ 2.5 it breaks down. We do not 
know at this stage the meaning of this critical value of the coupling constant, neither for 
our unphysical example of the scalar cubic interaction nor for other theories (e.g. $^). We 
even do not know whether the SDE indeed does not have a (formal) solution for larger 
A or whether it is just a failure of our numerical method. We would only observe that 
a very similar behavior was found in 0, where a similar coupling was considered. 
In our case it is also possible (in a slightly modified formulation) to calculate the field 
renormalization function Z which is finite. It appears that very close to \crit (within 5 
per cent) this function goes through zero (another sign of pathological behavior). We 
are aware of the fact that the model is in fact not defined at all |Q, hence its study 
should be viewed as purely methodological and one has to be very cautious in making 
any generalization. It might be nevertheless interesting that when one includes one loop 
corrections to the vertex function, Acrit drops by a factor of about 2. 

Figure 1 shows the weight functions p{a) from below for A = 0.1,0.5, 1.0,2.2,2.5. The 
last line already shows wiggles, precursors of the numerical breakdown. Figure 2 displays 
the self-energies below threshold s = Am? from above for A = 0.25,1.0,1.5,2.2,10.0, 
compared with those obtained by Dyson (bubble) summations (for the last value A = 10.0 
we do not get a solution of the DSE). For small A's the Dyson summation approximates 
our full solutions rather well, for A close to \crit they deviate by up to 25 per cent. 

4. CONCLUSION 

We have developed and tested in a cubic scalar toy model the method of solving DSE's 
in Minkowski space. Several applications of the method are in progress: studies of the 
scalar $^ theory, of 3+1 QED in a quenched rainbow approximation, and of a theory with 
'ijj^^ip^ fermion-pseudoscalar coupling. The solution of SDEs in the described framework 
seems to be feasible for all of these models, even with more complicated structure of 
equations employed. When one tries to include more loops, much more painful algebraic 
manipulations has to be done to bring equations for weight functions into manageable 
form. We hope to prove by this further tests that this method is competitive, compared 
to solutions in Euclidean space. 
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Fig. 2 Curves in figures are described in the text. 
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